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ABSTRACT 



Aims. We investigate the one-zone SSC model of TeV blazars in the presence of electron acceleration. In this picture electrons reach 
a maximum energy where acceleration saturates from a combination of synchrotron and inverse Compton scattering losses. 
Methods. We solve the spatially averaged kinetic equations which describe the simultaneous evolution of particles and photons, 
obtaining the multi-wavelength spectrum as a function of time. 

Results. We apply the model to the rapid flare of Mrk 501 of July 9, 2005 as this was observed by the MAGIC telescope and obtain 
the relevant parameters for the pre-flare quasi steady state and the ones during the flare. We show that a hard lag flare can be obtained 
with parameters which lie well within the range already accepted for this source. Especially the choice of a high value of the Doppler 
factor seems to be necessary. 
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1. Introduction 

Blazars, a class of Active Galactic Nuclei, have multi- 
wavelength spectra showing two components: the first, which 
is generally believed to be produced from electron synchrotron 
radiation, peaks, according to the object, from the optical up to 
X-ray regime; the second, which probably is produced from in- 
verse Compton scattering of t he same relativistic ele ctrons either 
on the synchrotron photons dMaraschi et al.lll992) or on some 
other, external to the so urce, photon population jDermer et al.l 
ll992l;ISikoraetal.lll994 '). peaks in the GeV to TeV y-ray band. 
This emission is thought to be produced within a relativistic flow 
pointing at us and therefore to be Doppler boosted. 

Flaring activity is probably the major characteristic of these 
objects. This is especially true for the X-ray and TeV regimes 
which are produced by radiation of electrons close to their maxi- 
mum energy. Thus in the X-ray regime flares have b een observed 
on timescales of hours to days (iKataoka et al.ll2001[) while in the 



TeV regime the variability can be even faster (iGaidos et al.l 19961: 
[Albert et al. 2007; Aharonian et al. 2007). 

Recendy it was reported (lAlbert et al.ll2007l) fliat Mrk 501 
showed an episode of flaring activity in the TeV regime where 
hard y-rays lagged the soft ones by approximately four minutes. 
This flare combines therefore fast variability with spectral evo- 
lution. Flares showing a hard lag h ave been observed befo re in 
the X-ray regi me (see, for exam ple, iBrinkmann et al.ll2005l and, 
more recently, Sato et al. '2008"), in addition to flares showing a 
soft lag (Takahashi et al. 1996). However, such detailed spectral 
information was not possible to obtain so far in the y-rays. 

In general, soft lags can be interpreted as due to electron 
cooling. If electrons are either accelerated very fast to their max- 
imum energy or are injected as secondaries, then as the more en- 
ergetic electrons cool faster, they will produce first a variation 
in the hard photons followed by a variation in the softer ones. 



exhibiting thus a soft lag ("Kirk et al.' '1998' - hereafter KRM, 
[Kii-k & Mastichiadis 1999. Kusunose et al. .2000) . 

Hard lags, by necessity, require some sort of p article accel- 
eration - for a different altogether explanation see l AIbert et"al] 
(2008). If the timescale associated with it is faster than the cool- 
ing one, then the radiation emitted from particles accelerating 
from low to high energies shows a hard lag. However, as the 
population of electrons moving to high energies does not have 
enough time to cool before reaching the maximum energy, the 
amount of produced radiation is small and therefore the mecha- 
nism is rather inefficient. Thus hard lags will be more apparent 
close to the spectral cutoff, where, by definition, the radiative 
loss timescale is comparable to the acceleration one (KRM). 

In the present paper we examine the unique properties of 
the July 2005 flare of Mrk 501 within the context of the one- 
zone SSC model which has been modified as to include electron 
acceleration. Our goal is to determine whether the model can 
give reasonable fitting parameters for both the steady state and 
the flaring activity. In Sect.|2]we outline the model, in Sect. [3] we 
show the procedure of fitting both the quasi steady state and the 
flare and we conclude in Sect. ID 

2. The Model 

In order to model the emission of Mrk 501 we use a similar ap- 
proach to Mastichiadis & Kirk (1997, hereafter MK97), i.e. we 
solve simultaneously two time-dependent kinetic equations for 
the distribution functions of electrons and photons. The electron 
equation reads 
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The operators £. denote losses and escape from the system while 
Q denote injection and source terms. The unknown functions 
He are Uy are the differential number densities of electrons and 
photons respectively and the physical processes which are in- 
cluded in the kinetic equations are: (1) electron synchrotron ra- 
diation and synchrotron self absorption (denoted by the super- 
scripts "syn" and "ssa" respectively); (2) inverse Compton scat- 
tering (both in the Thomson and Klein-Nishina regimes) ("ics"); 
(3) photon-photon pair production ("yy")- 

The physic al processes listed above are taken into account 
as described in iMastichiadis & KirkI (Il995h (hereafter MK95), 
where a detailed description of the functional form of the above 
operators is given. The current version of the code uses also 
some improvements over MK95. Thus for synchrotron radia- 
tion relativis tic electrons emit the full photon spectrum (see, 
e.g. Blument hal & Gouldl 19701) instead of the delta-function ap- 
proximation used in MK95. Also, for inverse Compton scat- 
tering, while the electron cooling still uses the technique de- 
scribed in MK95, the electro n emissivity uses relation (2.48) of 
iBlumenthal & Gouldl (Il970l) . Numerical tests have shown that 
this approach balances electron energy losses and total photon 
radiated power to within 90% of each other. 

Furthermore, departing from the approach of MK97, we im- 
plement an acceleration term in the electron kinetic equation 
which is characterized by an appropriate timescale (face) and 
is accompanied by a term which describes particle injection at 
some low energy (first term in RHS of Eq. [T] where Q{t) is the 
rate of electrons which are injected at low energies yQin^c^). This 
modification allows us to follow particles as they accelerate from 
low to high energies. Note that this approach is similar to the one 
taken in KRM. However, the present work differs from KRM in 
that we adopt here a one-zone model. 

There are six parameters that are required to specify the 
source in a stationary state. These include 

1. The Doppler factor (5 = [r(l -yScosfl)]-'. 

2. The radius R of the source (or, equivalently, the crossing 
time in the rest frame of the source - Rjc). 

3. The magnetic field strength B. 

4. The acceleration timescale face, which by the setup of the 
problem must obey the relation face ^ fcr- 

5. The timescale of particle escape of the system fesc- 

6. The rate of injected electrons - we note, however, that 
the solution turns to be largely independent of the exact choice 
of yo as long as this is not larger than 10. 

As was shown in KRM this prescription (under the assump- 
tion of synchrotron radiation losses only) leads to an electron 
distribution function which in steady state reads 

(3) 

r Tmax / 

for yo < y < ymax, where ymax is the Lorentz factor at which 
electron energy losses balance acceleration. For example, in the 
pure synchrotron case, 

Tmax = (jSsfacc)"' (4) 

where = I'^rCg^j^p- with crj the Thomson cross section. 
However, the present approach incorporates, in addition to syn- 
chrotron, SSC losses which render the derivation of an analytic 
solution impossible due to complications arising from the Klein- 
Nishina limit. One further notes that for face and fesc both in- 
dependent of energy, as it was assumed in deriving the above 
solution, the electron distribution function is a power law of in- 
dex s = -2 - (face - fese)/resc (as long as y « ymax). We note 



also that tests on the code in the pure synchrotron loss case have 
shown that the electron distribution above y,^ax does not drop to 
zero abruptly, but rather, due to numerical diffusion, a very steep 
power law is produced. 

The above description changes during a flare: Assume that 
the system has reached some stationary state. If this is per- 
turbed in some way, i.e. by injecting an increased number of 
particles in the acceleration mechanism for some time inter- 
val At, then a wave of fresh particles will move to high ener- 
gies. Assuming that the episode starts at some instant fo, then 
at each time f > fo + Ar the fresh particles will have Lorentz 

factors yfl,min(0 < y < yfl,max(0 with yfl,,nax(f) = yoe^'-'o^/'"'^' 

and yfl,mm(0 = -yQe('-'o-AT)//„cc 'pjjjg relation holds as long as 
yfl.max ^ ymax- As the time evolving particle distribution will 
have a higher amplitude than the steady state one, this will cause 
a flare in photons which will relax back to the pre-flare state once 
particles of Lorentz factor yfi,min start becoming of order ymax- It 
is interesting to note that as losses do not come solely from syn- 
chrotron radiation but from SSC as well, the flaring event could 
have an impact on the determination of ymax- More specifically, 
the presence of extra photons during a flare in the system will 
increase the total electron energy losses and, depending on the 
specific conditions, could cause ymax to drop. 

3. The July 9 2005 flare 

We begin by using the model described in the previous sect ion to 
fit the X/TeV data as given in Fig. 21 of I Albert et all (l2007i) (July 
9, 2005, observations). We have solved the set of stiff" diff'erential 
equations ([TJ and (|2]i using the numerical techniques as these 
were described and tested in MK95. The parameters used for the 
steady state fit are 7? = 1.5 10'"^ cm, 6 - 6Q, B - 0.5 G, face = 
3fci-, fesc = 4.17fcr, yo = 10°°^ and e = 5 10^ cm-^sec-'erg-'. 
The cosmological parameters used are Hq - 70 km s"' Mpc"', 
Qa = 0.7 and Q = 0.3. The spectrum is shown with solid line 
in Fig.[T] The resulting electron distribution function is a power- 
law of slope s - -1 .7 up to an energy ymax =2 10^. 

Perturbing the steady state as given above we found that a 
change in Q{t) always produces a hard lag flare as the one ob- 
served. However, this method of simulating a flaring activity 
causes the flux in each ene rgy band to increase approxima tely 
by the same amplitude - see IMastichiadis & MoraitisI (|2009|) . If, 
as the observations seem to suggest, the flare is becoming harder 
as well, i.e. more flux is emitted in the higher energy bands, then 
in order to get a fit to the TeV lightcurve one needs, in addition, 
to decrease face and/or B during the flaring episode. This can be 
explained from an inspection of rel. (4). A spectral hardening of 
the flare requires an increase of ymax during the episode and this 
can be achieved, within the context of the present model, only 
by reducing one, or both, of the aforementioned parameters. 

Figures 2 and 3 depict the lightcurves resulting from such 
a flaring episode. Here the observed flare comes from an im- 
pulsive change of Q by a factor of ~ 13 for Ar = Ifd and a 
decrease of face and Z? by a factor of 1 .7 for At - 30fcr which is 
approximately equal to the duration of the flaring episode, i.e. to 
the time needed for the injected particles to reach ymax- In order 
to keep the spectral slope unchanged (see rel. 3) we also vary 
fesc by the same factor Fig. 2 shows the lightcurve at the lowest 
energy range (0.15-0.25 TeV) of the MAGIC telescope, while 
Fig. 3 shows the corresponding lightcurve at the highest energy 
range (1.2-10 TeV). The time lag between the two energy bands 
is about three minutes. 

The chosen flaring episode is too short to approach the 
stationary state corresponding to the new increased injection. 
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Fig. 1. Fit of the multi-wavelength spectrum of Mrk 501 during 
the flare discussed in the text. Solid curve i s a fit to an assume d 
steady (or pre-flare) state as discussed by [Albert et"an (|2007|) . 
Parameters are as given in the text. The long dashed curve shows 
the spectrum at the instance where the fluence is maximum at 
the hard TeV band. Th e TeV spectra were de-absorbed as in 
iKonopeUco et al. I (l200l . 
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Fig. 2. Time evolution of the flare described in the text in the low 
(0.15-0.25 TeV) TeV band. Time is as measured by an observer 
in the lab frame. Observations were taken from lAlbert et alJ 
(I2007h . 

Fig-fflshows, in addition to the pre-flare steady state, the multi- 
wavelength spectra at the instance where the fluence is maxi- 
mum at the hard TeV band (long dashed line). 

Fig. |4] shows the the derived TeV hardness ratio (defined as 
F{1 - 10TeV)/F(0.25 - ITeV)) versus the TeV flux. This, as 
expected, is anti-clockwise. 

4. Summary/Discussion 

In the present paper we have attempted to explain the Mrk 501 
rapid variability as observed by the MAGIC telescope in 2005 
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Fig. 3. Time evolution of the flare described in the text in the high 
(1.2-10 TeV) TeV band. Time is as measur ed by an observer i n 
the lab frame. Observations were taken from i Albert et al.l ( l2007h . 
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Fig. 4. Hardness ratio F(l - 10TeV)/F(0.25 - ITeV) versus flux 
F{> 0.25TeV) for the flare described in the text. The curve ro- 
tates with time in an anti-clockwise manner. 



with a one-zone model which includes particle acceleration plus 
radiation. According to this particles are injected into some ac- 
celeration mechanism at low energies and subsequently they ac- 
celerated to high energies. They also radiate synchrotron and 
SSC radiation, thus losing energy. The acceleration process sat- 
urates eventually at some particle energy where it is balanced by 
the radiative losses. In this picture, the basic mechanism f or pro- 
ducin g hard lags as the one observed by MAGIC ( Alber tet al.l 
12007b . could be a brief episode of enhanced particle injection. 
As the freshly injected particles move to high energies they will 
radiate first soft photons and later harder ones, thus creating a 
hard lag flare. This episode occurs on top of a quasi steady state 
that the particles were assumed to have achieved in the pre-flare 
state and to which they return in post-flare. Thus according to 
this picture, MAGIC has caught electron acceleration at work. 
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While hard lags can be explained by increasing the parti- 
cle injection at low energies, the observed spectral hardening 
of the flare requires that, in addition, the acceleration timescale 
and/or the magnetic field strength should be reduced over their 
pre-flare values. Even if variation in one of the above parameters 
can produce a spectral hardening, varying both simultaneously 
improves the lightcurve fit. Given, however, the complexity of 
the overall flaring activity of TeV blazars, we consider the issue 
to be open, requiring more detailed observations. 

The parameters used h ere are quite similar to the ones used 
in iKonopelko et all (l2003h . Especially we note that in both fits 
there was a need for a high value of the Doppler factor 6. In 
the present case this c hoice is almost obliga tory given the fast 
variability of the flare dBegeknan et al 

The model presented here is an extension of the KRM model 
in the sense that it uses, in addition to synchrotron, SSC losses 
to balance acceleration. However, contrary to KRM, it is an one- 
zone model, thus it is implicitly assumed that particles which 
escape do not radiate. This can be justified only if the magnetic 
field strength in the region where the particles escape is much 
lower than in the radiating zone. Also SSC losses are expected 
to be reduced away from the radiating zo ne as the s ynchrotron 
photon energy density drops away from it (lGouldll979,) . Finally, 
we should mention that despite the fact that the electron equa- 
tion (Eq. [TJ is similar to the one used in the kinetic description 
of diffusive shock acceleration (e.g. Ball & Kirk 1992), we do 
not specify here a detailed method for particle acceleration. It 
is obvious though that simultaneous variations of more than one 
parameters, as is the case here, need to be interconnected with a 
specific theory for acceleration. 

A prediction of this model is that a flare due to enhanced 
particle injection causes a hard la g behavior in the X-ray regime 
as well (Mastichiadis & Moraitisi2009 ). However, in non-steady 
situations, as it is probably the case with flares, one expects to 
get complicated relations between the synchrotron and the SSC 
component (see Katarzvnski et al. 2005b. 

Recently Bednarek & Wagner ( 2008) have proposed another 
explanation for the hard lag by allowing the radiating blob to ac- 
celerate during the flare. The authors predict that, in such a case, 
the delay between 20 and 200 GeV should be approximately one 
order of magnitude longer than the observed delay in the TeV 
band. On the other hand, our present model predicts that such 
fast flares should have the same time lag across the y-ray band. 
Thus coordinated future multi-wavelength observations involv- 
ing GLAST and Cherenkov detectors could disentangle whether 
the lag is due to acceleration of the blob itself or whether it is due 
to a brief episode of particle acceleration within a blob which 
moves at constant Lorentz factor T. 
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